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Quantum fluctuations of a nonminimally coupled scalar field in D-dimensional homogeneous and 
isotropic background are calculated within the operator formalism in curved models with time 
evolutions of the scale factor that allow smooth transitions between contracting and expanding and 
between decelerating and accelerating regimes. The coincident propagator is derived and used to 
compute the one-loop backreaction from the scalar field. The inflationary infrared divergences are 
absent in Bunch-Davies vacuum when taking into account a preceding cosmological era or spatial 
curvature which can be either positive or negative. It is found that asymptotically, the backreaction 
energy density in the minimally coupled case grows logarithmically with the scale factor in quasi-de 
Sitter space, and in a class of models decays in slow-roll inflation and grows as a power-law during 
^ ■ super-inflation. The backreaction increases generically in a contracting phase or in the presence of 

| a negative nonminimal coupling. The effects of the coupling and renormalization scale upon the 

quantum fluctuations together with the novel features due to nontrivial time evolution and spatial 
curvature are clarified with exact solutions and numerical examples. 
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O , 1 I. INTRODUCTION 
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At the largest scales, the universe seems very nearly homogeneous and isotropic. Thus its evolution can be described 
by a single scale factor a(rj), whose time evolution is parameterized by the conformal time rj. The metric, originally 
studied by Friedmann, Lemaitre, Robertson and Walker (FLRW), can be written as 



i , > ( , o dr 



g^dx»dx v = a 2 ^) [ -d V 2 + 1Kr2 + r 2 d6 2 + r 2 sin 2 6d(j) 2 ) . (1) 



The curvature of the spatial hypersurfaces is given by the constant K . The current cosmological measurements [l| give 
a very large lower limit on the radius R c — \K\~? of the curvature of the universe, R c > 22ft, -1 Gpc (R c > 33ft -1 Gpc) 
for a universe with positively (negatively) curved spatial sections, where ft = 0.705 ± 0.013 is the Hubble parameter 
evaluated today in units of 100 km/s/Mpc. In the present work we will work in terms of the conformal Hubble rate 
defined as 



Hfa) = ^log(ofa)). (2) 

The leading paradigm explaining the observed approximate isotropy, homogeneity and flatness of the universe is 
inflation Q . Remarkably, an inflationary theory provides also an origin for the observed structures in the universe as 
a result of quantum fluctuations in expanding space 



Scalar fields have been ubiquitous in the studies of this phenomenon. The behavior of a massless scalar field 
with a possible non-minimal coupling to the Ricci scalar on expanding spaces Q has been studied extensively. This is 
relevant also for understanding of higher spin fields. In particular the graviton can be related to the kinetic operator 
of a scalar field. Apart from some tensorial structure, the graviton propagator can be written in terms of propagators 
of massless scalar fields 0] also in quasi-de Sitter space Q , and similar observations can be made for vector fields 
[?J or antisymmetric tensor fields 0. It is also shown that a three-form inflation in some limits allows an equivalent 
scalar field description 0, [l(| ■ Therefore the understanding of the massless scalar field can be translated to more 
general fields as well. 
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The key property of quantum behavior of fields in an expanding spacetime is the particle production that stems 
from the presence of a horizon. The particles are created mostly in the infrared, and it is this effect that leads 
for example to the creation of primordial density fluctuations. This also raises the issue of the backreaction of the 
fluctuations on the background spacetime [ll], [l2j]. When one chooses the vacuum to contain only purely positive 
frequency modes (the so called Bunch-Davies vacuum [l3[ ), the expectation value of the two point correlator for this 
vacuum diverges in the infrared. Unlike the ultraviolet divergences, which can be removed by standard procedures, 
these infrared divergences show that the system under study is pathological. Hence, it is clear that the Bunch-Davies 
vacuurrQ cannot describe a physically sensible state of an exponentially expanding spacetime. Possibilities that 
emerge are to consider an alternative vacuum state with less singular superhorizon modes [l5| or a compact spatial 
manifold where those modes are initially absent (l6j . 

Many studies though concern the special limit of exactly exponential inflation, i.e. de Sitter background fl7HT9T ] . It is 
natural to consider more general evolutions to understand how the predictions from more realistic assumptions may 
differ from the maximally symmetric case. Recently the scalar and graviton propagators have been derived for flat 
FLRW universe with constant deceleration parameter [20l | . The infrared divergencies persist for a large and physically 
relevant range of the deceleration parameter, in particular for accelerating regime if the field is minimally coupled. 
However, in Ref. [21| it was shown that the infrared can be regulated by taking into account a decelerating epoch 
preceding inflation. The perturbation modes were matched at a transition from radiation domination to inflation, 
and the results agreed qualitatively with the ones obtained by using a cut-off procedure in Ref. [22j , when the effects 
of the particle creation due to the sudden transition approximation were properly identified. 

In this paper we consider universes with a smooth transition to an inflationary period. This is achieved by con- 
sidering more general background evolution specified by three parameters, which still allow analytic expressions for 
the solutions of the mode functions. Furthermore, we also include the possibility of nonzero spatial curvature. It is 
natural to consider that the universe was not exactly flat but the curvature radius was inflated beyond the present 
observational limits. In such a case, we cannot observe the curvature directly but we obtain consistent results from 
the calculation of fluctuations. In a similar way, by taking into account a bouncing or a decelerating phase preceding 
inflation will rather generically eliminate the infrared infinities in the quantum fluctuations. Thus it turns out that 
in the more complete models featuring either some curvature or realistic time-evolution, the inflationary spectra are 
regular. As some details of the spectra depend on the physical means of regularization, there is also a hope to obtain 
some information on pre-inflationary cosmology. 

Having then well-defined and regular quantum fluctuations allows to investigate their physical effects. It would 
be of interest to extend the study to the graviton. In this work we however focus on the leading order corrections 
from the massless scalar field in particular during inflation. For this purpose we derive the one-loop expectation 
value of the induced stress energy tensor for a scalar field. Our results confirm the well-known leading logarithm 
behavior in quasi-de Sitter inflation 0, 0, HH (see the interesting recent discussions on the physical relevance of such 
logarithms and their implementation within the 5N formalism [24], H3> a l so [Hi H3 an d the review [28[ .) If the 
backreaction grows, the perturbative treatment will eventually become invalid. There are suggestions that such may 
signal analogous breakdown of locality in quantum gravity in cosmological as in black hole physics [29]. We find 
that the backreaction generically increases also when the scalar field has a negative coupling to the Ricci curvature, 
in agreement with the conclusion reached in Ref. 21], or when the universe is contracting or super-accelerating. In 
most cases the sign of the induced energy density can change dynamically and depends on the parameters of the model. 

We set up our notation with some preliminary considerations in section [TT] where we also describe the background. 
The propagator is derived in section HTT1 using the operator formalism. For the above mentioned reasons, particular 
care is taken to obtain physically regular infrared, since the divergences (except logarithmic ones) would be absorbed 
by dimensional regularization. The ultraviolet is renormalized by counterterms. In section IIVI we then apply the 
propagator to compute the backreaction from quantum fluctuations. In particular, the asymptotically dominating 
contribution at inflation is extracted. Finally, we briefly discuss the conclusions in section [Vj Some details of the 
D-dimensional curved FLRW models are confined to the appendices: in appendix [Al we construct the curvature from 
the metric and derive the field equations, in appendix [B] we state the spherical Laplacian operator and discuss its 
eigenmodes. Appendix [Cl contains explicit expressions for the propagator in a few special cases. 



Here our definition of the Bunch-Davies vacuum differs from the one used in e.g. Ref. [14)1 , where the positive frequency solutions are 
defined at the limit rj — > — oo. For calculational purposes it is impractical to take this limit when the wavenumbers k considered range 
from to oo, and we pick the positive frequency modes at a finite rjg. 
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II. SCALAR FIELD IN EXPANDING SPACE 

Here, in III Al we will sketch the usual quantization of a scalar field in FLRW universe, and propose the ansatz (fT5|) 
for its background evolution. We then describe this class of background evolutions in III Bl 

A. Scalar field 

Consider a massless non-minimally coupled scalar field $ in a D dimensions, 

S=\ f d D x^-g®{U-ZR)$. (3) 



2 . 

Here R is the Ricci scalar and g the determinant of the metric. The d'Alembertian operator is denoted as □ = V^V^. 
The field obeys a simple equation of motion 

(□-£R)$ = 0. (4) 

The stress energy tensor becomes 

2 S S 1 

T^ = —— = - -g^f + £ (GV - V^V, + 9lu ,U) $ 2 , (5) 

with the trace 

T = gTTp, = ((d$) 2 + dm 2 ) + (D- 1)£D$ 2 . (6) 

We shall be interested in deriving the Feynman propagator iA(x;x'), which satisfies the equation 

y /=g(a-£R)iA(x;x')=i5(x-x'). (7) 

For this purpose we work in the operator formalism, and eventually sum over the mode functions (solutions of the 
equation (j4]) in Fourier space) to compute the propagator. This way we obtain the complete answer, whereas a 
solution to the Eq. (J7J is unique only up to the homogeneous modes. 

The metric, regarded as a background, was introduced in Eq. ([1]). Due to isotropy, the Fourier components of 
the field can then be expanded with the creation and annihilation operators 6k and fe^ as 

<S> k (v)=Mv)bK + r k (v)bl k . (8) 
The convention we use for the transformation is 

F(x) = * / d D - l kF^, x) . (9) 

(27Tj 2 J 

In a closed universe the integral is more properly replaced by a sum. In flat universe the mode functions are simply 
plane waves, Q(k,x) = e , and some details of the curved cases are given in appendix 151 So the comoving 
wavelengths of perturbations are ^ a/k. The canonical commutation relations for the operators and the field, 



(2^) IJ - 1 ^- 1 (k-k), (10) 



$ k fa), a( V ) D - 2 ^l)'( V )\ = i{2^r-^-^ - k) , (11) 
imply that the Wronskian of the mode functions is given by 

MvWkYiv) - PkhWM = ™ 2 ~ D W ■ (12) 

Above and in the following, a prime indicates a derivative with respect to the conformal time, ' = d/ drj. From now 
on we will occasionally drop the subscripts k and arguments of r/. 

To present Eq. Q in the metric ([T]), we need the Ricci scalar (|A6I) from appendix [S] and the dAlembertian 
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operator (|B18|) from appendix [Bj Notably, the Laplacian part of the latter, when acting upon the Fourier modes, 
gives a combination reducing to (k 2 — K) in D = 4. By the end of the day, the mode functions ipk are found to satisfy 
the evolution equation 

d 2 +k 2 + (~(D - 2) + (D - (D - 2)K + M 2 (H) 

where we have lumped the Hubble-dependent terms into the definition 

M 2 {H) = -i (D - 2) (n' + ^ (D - 2) n 2 ^j +(D-1) (2U' + {D - 2) H 2 ) £ . (14) 

Because this allows us an analytic solution for the mode functions, we assume that 

Af 2 (H)=(V + i^) . (15) 

where y 2 and v are D-independent constants. This parameterization captures a wide range of interesting geometries 
in D = 4 as described in the following subsection. The generalization of the geometries to higher dimensions is by 
no means unique. One could add arbitrary terms proportional to (D — 4) into equation (1151) . but we chose the trivial 
£>-dependence for simplicity. Such choices in an analytic extension affect only the finite constants that appear in 
the results. Note that after we have dictated (fT5j) . the rest of the terms in the evolution equation (fT5)l are uniquely 
determined by the action ([3]) . The remaining ambiguity will appear in the choice of the -D-dependence of the initial 
conditions for the solutions to this equation. This ambiguity is fixed by the natural prescription we adopt in section 
IIII A( namely that in any D the modes are initially in their Bunch-Davies state. 



o. 



(13) 



B. Background 




FIG. 1. Evolution of the slow-roll parameter e as a function of the conformal time 77 in units of |7|. Left: Solutions to Eq. fj 16p 
with positive j 2 . The curves start from the beginning of the universe. Right: Solutions to Eq. (|16|) with negative 7 2 . The 
curvature dominated phase can be extrapolated to r\ — s> —00. 



Let us look at which kinds of background evolution the ansatz (fl"5j) allows. In D = 4, the evolution equation for 
the scale factor becomes 



We can set £ = for simplicity. A nonminimal coupling corresponds then to redefining the parameters as 

^^ = VT^' 7 ^ = 7t^T (17) 
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when £ < min(|, 2v 2 /i). One sees that though the asymptotic behaviors at large and small 77 are the same as in a 
cosmology with curvature and an energy component with the constant equation of state w = (3 ± 2u)/(3 =p 6v), the 
geometry of (|15p is not precisely reproduced by the latter configuration but generally requires a dynamical w. In the 
constant-w case, Eq. (fl"5|) would be replaced by 



— = iC 



3, 
2* 1 



3to)|Jf|5" 



(18) 



where S(x) = sinhx for negative and S(x) — sinx for positive curvature. Since S(x) ~ x + . . ., our results should 
qualitatively describe also the particular case of spatial curvature and a constant-!/; fluid as the energy sources. In 
more detail, the evolutions we have that lead to inflation can be classified into three qualitatively different types as 
we describe below. We will later need the scalar curvature for these geometries, which is 



R 
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(1 - 6£)a 2 



7 2 + 



- (1 - 6£)K 



(19) 



when D = 4. 



1. j 2 > 0: From radiation domination to inflation 

In the explicit examples here, we assume that v is a half-integer v > i but this assumption is not crucial to our 
conclusions as explained later. Let us first consider the case of real 7 > 0. Then the exact solution to (|16|) can be 
expressed in terms of the Bessel functions 



0.(77) = (-ly-h^l [ Cl J v ( V j) + c 2 Y v { V i)] , (20) 

where c\ and C2 are real constants. In an accelerating universe, 77 approaches zero from below. At late times, from the 
asymptotics of the Bessel function, we see that 0,(77) — > \c2\2 IJ ~ 1 T(i')\'yri\ 2 _I/ , so the solutions describe the evolution 
infinitely far in the future where the scale factor grows without limit when v > i. This power-law evolution implies 
that the deceleration parameter will be asymptotically a constant, which is given by v as e = (3 ± 2i/)/(l ± 2v). In 
this paper we consider the upper sign choice since we are interested in inflation, but it should be kept in mind that 
for example v — 3/2 can correspond to matter dominated era as well as de Sitter spacetime. We use 



(21) 



v? 

corresponding to the usual definition of the slow-roll parameter. It is sometimes convenient to refer also to the 
deceleration parameter q — e — 1. If the scale factor decelerates with respect to the cosmic time, q is positive, and in 
the accelerating case q < 0. These are what we refer to as accelerating or decelerating in this paper (one notes that 
acceleration with respect to the conformal time a" > corresponds to rather e < 2 and q < 1). 

In the early times, the solution will however be oscillating, as immediately seen from Eq. (1161) when £ < 1/6. 
This means we cannot extend the evolution infinitely far in the past, but there is a beginning to the universe. 
Inspecting the behavior of the solutions, we see that generically the evolution in the physical region (positive scale 
factor) begins with the slow- roll parameter near e = 2, corresponding to radiation domination. Numerical examples 
are plotted in the left panel of figure [TJ These correspond to such initial conditions that c\ = and thus 

K=^ tl ^M. (22) 
2?7 M7?7) 

We note that a possible feature is a period of super-inflation (negative slow roll parameter) during the transition from 
e fa 2 to the constant- e phase given by the first term in (|2"2")l . 



2. 7 2 < 0: From curvature domination to inflation 

Let us then consider an imaginary 7 in equation (|15[) . The early behavior then becomes exponentially growing in 
77 instead of oscillating. Then e = 1, corresponding to a curvature-dominated universe. The exact solution can then 
be written in terms of modified Bessel functions, and we can extend the evolution infinitely far in the past. We have 



a(v) = V\rn\[cil^(hv\) + C2K^(\jri\)} , 



(23) 
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where c\ and C2 are again real parameters. We plot some example evolutions corresponding to the initial condition 
c\ = in the right panel of figure [TJ This corresponds to the conformal Hubble parameter 

*~ + l7l%#". (24) 

We chose the normalization in such a way that future asymptotic expansion is the same as in the previous case. It 
turns out we must include curvature in the case of imaginary 7 to consistently regulate the infrared. The required 
amount is \K\ > \j 2 \ if K < and K > \-f 2 \/3 if K > 0, as will be clarified later. 




FIG. 2. Bouncing universes described by Eq. (|23[) when j 2 < and C2 = 10ci. Lowering [--y | corresponds to stretching the 77-axis. 
Left: The Hubble rate H = H/a for several values of v. The evolution begins from a contracting phase H < 0, experiences a 
bounce at H = and enters an accelerating phase. In case v = 3/2 the latter becomes de Sitter and H approaches a constant 
Ho. Right: The effective equation of state of the background fluid for the v — 3/2 case. In closed universe, the fluid obeys 
the null energy condition. This is not the case in open or flat universe (and the latter example is not infrared regular). 



3. ci,C2 > 0: Bouncing from contraction to inflation 



A third qualitatively different case arises when we include allow also a nonzero c\ in either the solution l|20p or 
the solution (|23[) . For the sake of clarity, we consider only the latter case here. Then we find a contracting phase as 
rj — > —00. Thus the inflationary epoch was preceded by a cosmological bounce. Such scenarios have received lot of 
attention [3(|, because they may resolve the cosmological singularity and extend the evolution to sort of a pre-Big 
Bang phase. A bouncing universe can be geodesically complete unlike an eternally expanding one [3lj ] - 

Few examples are plotted in figure [2j In the left panel we show the Hubble rate H/a. The right panel depicts 
the effective equation of state wb for the energy component that is the source for the Hubble rate. It is well known 
that in general relativity, only null-energy condition violating energy sources can drive a bounce in a flat universe, 
whereas in the presence of curvature it can be sufficient to violate the strong energy condition. One may however 
obtain bounces in flat setting without introducing ghosts in generalized gravity theor ies, recent progress has been 
made within the first order formalism [32I |33T | and for string-inspired nonlocal gravity [34 . l35| . Using the equations 
in appendix El we have deduced the equation of state wb of the background fluid in the present case. The general 
expressions are rather cumbersome, but let us quote the result in the special case v = 3/2: 

cf [t (t ({t + 1) 2 k + 2) + 6) + 3] + 2cic [i 4 (2 + «)- t 2 {n - 4) - 3] e 2t + c 2 [t (t ((< - l) 2 n + 2) - 6) + 3] e 4 * 

WB ~ 3c 2 [t(t + l)(t(t + j> + 2) + 1] + 6cic [i 4 ( K - 2) - t 2 n - 1] e 2t + 3c 2 [(t - 1)£((£ - l)tn + 2) + 1] e 4t ' 

(25) 

where 

K 

*=|7|?7, K = l c = ci + 7tc 2 . (26) 
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Y>0 ■?=<) y 2 <0 




FIG. 3. The colored areas in the curvature-coupling plane are excluded because of the infrared divergence. The horizontal lines 
indicate the conformal coupling £ = |. Left: The cases 7 2 > (curvature in units of 7 2 ) whose background evolutions were 
discussed in MB ll The minimally coupled models £ = are always regular. Middle: The case 7 = 0, which corresponds to 
constant deceleration spacetimes. Note that the flat models K = are excluded. Right: The cases 7 2 < (curvature in units 
of — 7 2 ) whose background evolutions were discussed in III B ll and lIIB~3l The minimally coupled models £ = require K < j 2 
or^>|| 7 | 2 . 



The effective equation of state is plotted in the right panel of figure [3] In our case, it turns out that for general v 
and C1/C2, when K > 0, the background fluid responsible for the bounce indeed typically respects the null energy 
condition. Thus we could reproduce this behaviour with e.g. a canonical scalar field. Obviously, with negative 
curvature one has to violate the null energy condition since ps ~ T~L 2 + K. Except when K > 0, the collapsing phase 
must be driven by a negative energy source. In the simplest form this is a negative cosmological constant, which 
however has to go through a phase transition to recover an observationally allowed vacuum density at later epochs; 
such scenarios have also been elaborated recently [3(1 137[ . Here it is natural to consider the bounce to occur in a 
closed universe where the collapsing phase is driven by supercritical density of matter. To generate a bounce, this 
density needs to violate only the strong energy condition. In the following we find that the inflationary infrared 
divergences are absent in such models. 

To close this section, we note that when v is not a half-integer, the solutions for negative 7 2 are written then 
in terms of the modified Bessel functions, and the solutions for imaginary 7 2 in terms of Bessel functions. This is 
the opposite to the above cases in subsections III B II and III B 21 Furthermore, with a positive nonminimal coupling 
£ > 1/6, the solutions once again switch their roles. 



III. SCALAR PROPAGATOR FOR QUANTUM FLUCTUATIONS 



Here we first derive the scalar propagator and evaluate it at the coincidence limit in IIII Al The cases where this is 
due to infinities lost in the analytic continuation are identified and discarded (see figure [3]) . We then renormalize the 
ultraviolet divergencies in a standard way by adding suitable counterterms in IIII Bl In section IPl! CI we show that the 
previous results are recovered in the appropriate limits. 



A. Infrared-regular sum over the modes 

The solution to the mode equation (|13[) with the ansatz (|15|) can be written in terms of the Hankel functions, 



Mv) = a 2( ,-n)\Mr 



TT\T)\ 



(27) 



where we have defined 



2 2 
Ik = 1 



~(D- 2) -(£>-!)£ 



(D-2)K, 



(28) 



which plays the role of a shifted j 2 in the presence of curvature. While j 2 specifies the evolution of the background, 
it is the 7^ which determines the properties of the fluctuations. The normalization from the Wronskian (|12[) is such 
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that \a\ 2 — |/3| 2 = 1. The Bunch-Davies vacuum corresponds to the choice a = 1, /3 = 0, which we make here. The 
time-ordered Feynman propagator, obeying Eq. ([7]), for the vacuum state |f2) is 

iA(x; x') = (Q\e(r) - t/)$(x)${x') + 9(rf - r?)$(a/)$(x)|fi) 
dP^k 



(27T) 



D- 



jQ^k^Qik^ieiv-^ivWi^ + oiv' -vmv'Wiv)} ■ (29) 



In the second line we recalled the expansion © and used the definition of the vacuum &k|fi) = 0. In the case of 
closed universe, one should, more rigorously, replace the integral by a sum over the modes as mentioned in appendix |B| 

We shall need the propagator at the coincident limit. At coincidence the mode functions Q(k, x) reduce to plane 
waves. After performing the angular integral in (1291) and inserting the solution (|27[) we get 

l A(x;x) = f°° dkk D - 2 HW(Jk* + i> K \r,\)H<?>(Jk* + , (30) 

2 L, Tr—2-r(±Y j -) Jc v 



where C is the cut-off that equals DVK/2 when K > and vanishes C — when K < 0. At this point we assume 
that v is a half-integer. Then we may employ the finite series presentation ^ v 2 

V 7 

and obtain 



V 7tz *— ' n!l (z/ + £ — n) 



, , , ° 2 -° r dfcfc ^- 2 v <m "" r(,/ + * + n)I> + * + m) (fc2 + 7p-^ (l+ " +m) f321 

1 ' ' ~ 2°-^T( V) Jc n j^ Q n\m>nv+\~n)T{v + \- m ) (2|t,|)«+™ ' ^ 

The benefit is that the integral can be performed analytically. When the radius of convergence is nonzero, it can be 
shown that the result applies well for general v [2l|, however in the present case the analytic extension seems more 
difficult. From the form of the integral we see that divergence is expected at the upper limit when D — 2 > n + to, 
and that terms in the sum with odd n + to vanish. In fact the double sum can be reorganized into a single sum as 



roc _2 F ( i + n)T{ 2 v + 1 + 2n)T ^ + i _ ra ) ( fc 2 + 7 |)~|- 

2£>-i 7r f r(^I) Jc ^ n!r(2^+ 1 - 2n)T(v + 1 + n) (4t7) 2 ™ 



n 2-D roo - 2 



= « 2 - D f°° B _ 2 v 5 r(i + n)r(, + i + n) (fc 2 + 7 j)-i-" 

2 c-i vr fr(V)^ ^ n!I>+|-n) r? 2 " 

In the second line we used the Legendre duplication formula. We note that j%- > —C 2 is needed for the regularity 
at the lower limit of the integral. However, the results turn out to be analytically extendable to arbitrary imaginary 
jK, but this is due to the automatic subtraction at play in the dimensional regularizatior@. This is why we discard 
the 7^- < — C 2 cases and can restrict to real jk in flat and open universes and consider negative 7^ possible only in 
closed universes. The allowed spacetimes we are left with are then the following for the minimally coupled case: 

• 7 2 > 0: This case is always regular. 

• 7 2 = 0: Spacetimes with constant deceleration are regular when the spatial sections are curved, X^0. 

• 7 2 < 0: This case is regular for open universes with |ijT| > mz^p or closed universes with K > -jyTj- 

These are the cases in which the integral in Q33[) yields finite results without hiding the infrared power-law divergences 
by dimensional regularization. The regularity conditions in the presence of a nonminimal coupling are illustrated in 
figure |3l Here we see that it is the combination of the background evolution, given by 7, and the curvature, given by 



2 The 1/2-terms in the Gamma-factorials have a different sign in Ref. [38| 4th edition in Eq. (8.466), but Eq. (8.451) is correct. 

3 The incorrect use of dimensional regularization has been noted to hide power-law divergences in constant-e spacetimes \'2'M and the 
breaking of de Sitter invariance of propagators in e = spacetimes [39t . 
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K, which determines the physically reasonable cases. We also note that both positive or negative spatial curvature 
can regulate the infrared, but by a different mechanisms. Negative K contributes to the evolution of the fluctuation 
modes, mimicking a real 7, and can thus render these modes less singular at very large scales. On the other hand, 
with a positive K the very large scale modes are simply absent, because they cannot be excited in a closed universe. 



Keeping in mind the restriction to the above cases, we perform the integral in (|33|) to obtain 
iA(x; x) 



C D - 2 T(\ + n)T{v + \+n) 



^-2 2 o-i 7r #r(^i) ^ n!(2(l + n) - D)T(v + \-n) 

( D 1 D 7r-\ 1 

X ^{ 1 -2 +n '2 + ^ 2 -2 +n '-^)jdi^' (34) 

in terms of the hypergeometric functions with arguments proportional to square of the inverse cut-off. The result 
diverges in the limit jk — > if v > |. Thus by introducing 7^ we have indeed regularized the infrared. However, an 
ultraviolet divergence remains. We will deal with that in the following subsection. 

First we note that it is useful to rewrite this result as 

1k~ 2 r (! - f + + \ + ") 1 



iA(x; x) 



22^* ^ nW(p + \ - n) {lKV? n 



IK J r(i-f +n)r(^±i) ' *V2 V "2 ' '2^ " 7^ 



(35) 



The first line is the limiting value of the above (|3~4")) as C — > 0, and is thus the complete result for flat and open 
universes. The second line deviates from unity in the presence of a nonzero cut-off C, and in particular it takes into 
account the absence of modes with k < in closed universes. Especially in the case of vanishing cut-off this 

expression is more practical than (|34[) . However, we should stress that the Gauss hypergeometric functions appearing 
in (|35j) have a branch cut at imaginary jk and in that case the use of the expression (|35j) requires some extra care. The 
expression ()34j) is straightforwardly applicable in the regime jj^/C > — 1 covering precisely the physically regulated 
cases. The expressions are equivalent, but for practical purposes (|34[) should be used when C 2 3> 7^, and ([33)) when 
C 2 <C 7^ (and in particular when C = 0). 



B. Regularizing the ultraviolet by counter-terms 



Let us first look at the case ([35| . The first two nonvanishing terms in the series (|35|) contain the divergent piece: 



,2-15 



iA uv (x;x) 



D-4 



IK 



r(i-f)[8 7 ^ 2 + (^-2)(i-4, 2 )] 



1 1 



r(i(£> + i)) 



8-yW x 2F X ( -, ^(D - 1); \{D + 1); -gl 
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7| 



t^'' J I) - ^1 ( |,i(i?-l);i(X> + l);-^- 



(36) 



where /1 is an arbitrary renormalization scale with the dimension of mass. The large-scale cut-off is not relevant for 
the ultraviolet divergence, and is thus contained C — part of the above propagator. This is the first line of (|36|) . 
The hypergeometric functions in the two following lines are regular in D = 4 and can be straightforwardly evaluated 
there. This becomes, near D — 4, 



S- D iA uv (x;x) 



1 - Av 



2 r 



1 - 4v 2 + 4(7k??) 2 

327r 2 aV(L>-4) ^ 64AV 



7b + log 
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1 -^- l0g (iJ^) +2 7| (1 - 5e) 
C(l~4^ 2 + 4( 7 2 f + C 2 )r; 2 ) 
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IK 



0{D - 4) .(37) 
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We have multiplied the propagator by /x 4_£) to obtain the correct mass dimension off D = 4. Here jk is evaluated in 
D = 4. We denote the Euler-Mascheroni constant by 7^ and ip(x) — 4- log (r(z)) is the digamma function. The last 
line in (|37|) contains the contribution from the hypergeometric functions that vanishes when C = 0. In the minimal 
subtraction scheme, all but the first term in Eq. Q37[) contribute to the renormalized propagator. One may write an 
equivalent expression either by directly expanding the first two terms in the series in ([34]) about D — 4 or by using 
the previous result ([37]) . We obtain 



iA uv {x;x) 
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1 - Av 2 + 4( 7J f t?) 2 
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8C 2 



C 2 



(38) 



up to terms vanishing in D = 4. This form is more convenient to use when C is large compared to \^k\ 



We will be interested in computing the one-loop expectation value for the trace of the scalar field energy mo- 
mentum tensor ([5]). The value of the trace © is given by, after using the Eq. ([7]) and the equation of motion 

(□-eR)o> = o, 



(0|T|0) 



D 



{D-l)i )DiA(x;x) 



(39) 



At this point we note that the ultraviolet divergence (|3"T)) is proportional to the scalar curvature (TT91 . This is what we 
expect, since the one-loop correction is known to be renormalizable with simple geometric counterterms. Generally, 
we may consider the two quadratic geometric counterterms (as the Riemann squared in D — 4 is a linear combination 
of these) 



This lagrangian contributes to the trace as 

T ct = [-4(D - l)ai + a 2 D] OR+(D- 4) (aiR 2 
Thus the (infinite parts of the) coefficients that renormalize the theory are 



a 2 Ra V R^) 
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(40) 
(41) 

(42) 



As jk is related to infrared physics and doesn't enter into this expression, it coincides with the result in 21]. The 
total trace from the divergent part of the scalar field stress energy ((39)) and the counter terms (l4Tj) becomes 

(n\T div + r ct \n) = f 7 | + z^-) + 12(D ~ 4)Q2 [2h 4 + m 2 W + 2{n'f - k(ah 2 + m' - ak)} 



32n 2 a 
24 



if 



3a{ m - af " ) (H'" - m 2 W - 2(H') 2 + 2H'K 



Jin 



(43) 

The second line contains the contributions from the possible finite parts of the coefficients a.\ and a 2 - A simple choice 
for renormalization is to set 012 and the finite parts of both coefficients to zero. The latter corresponds to the minimal 
subtraction scheme in dimensional regularization. We will use this choice in the numerical examples that follow in 
section IIVI 



C. de Sitter and quasi-de Sitter limits 



1. de Sitter limit 



It is known for a long time that in a locally de Sitter space the propagator of a minimally coupled massless scalar 
at the coincidence limit contains a logarithmically growing part [3, 4, 23] . To check the consistency of our results with 



4 We note that in Il37|l an imaginary 7u inside the logarithm would require an imaginary fi for consistency. This would spoil the unitarity 
of the counter term lagrangian 1 140 I I. which is one way to see that the 7^- < cases are unphysical in flat and open universes. However, 
we can allow an imaginary -yx m a closed universe. As argued above, the propagator (134 6 is then devoid of infrared singularities given 
C > \ik\, and when expanded about D = 4 in Eq. 1381 . is also perfectly compatible with a real renormalization scale fi. One should 
just note that in the case of imaginary 7x, the different branch-cut should be chosen in Eqs. H35I37II . 
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the previous literature, we consider the asymptotic limit when v = 3/2 and 



1 



(44) 



At small enough \r]\, i.e. after a sufficient inflationary expansion, this becomes an excellent approximation. Setting 
the cut-off to zero, the propagator (|35[) becomes, near D = 4, 



where the constant £ is defined as 
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(45) 
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1G 



(46) 



Eq. (|45|) agrees with e.g. the result derived in Ref. [40]. In the first line, the prefactor is the familiar square of 
the Hubble rate per 47r, then we have a divergence, and the secular growth which is logarithmic in the scale factor. 
The following constant depends on the renormalization scale and other parameters. In the case of nonzero cut-off, 
this constant would have a different form. The second line represents a correction which decays as 1/a 2 , and can 
be neglected at sufficiently late times (the constant coefficients would again be modified if there was a cut-off). All 
such redshifting terms can be absorbed into a suitable modification of the initial state [4l| . Our results consistently 
generalize the previous coincident de Sitter propagators into the case where the initial (Bunch-Davies) state does not 
have to be set in the exact de Sitter background. 



2. Quasi de Sitter limit 



We have essentially three parameters which determine the scale factor expansion: 7 and v specify the evolution 
equation (fTl)]) and the ratio c-ijc\ of the constants in the solutions discussed in section Hi Bl specify a particular 
solution to this equation. Qualitatively different classes of spacetimes emerge when the sign of 7 2 is changed or one 
of the coefficients c\ or C2 is switched on or off, as detailed in section III Bl However, we were able to carry out the 
computations analytically only for a discrete set of ^-values. This seems to restrict us to consider only a limited set of 
inflationary models. In particular, to study spacetimes which become asymptotically infinitesimally close to de Sitter 
would require to treat v as a continuous parameter. Analytical extension of the result (|34[) to arbitrary v though 
does not seem straightforward as the coefficients in the series diverge factorially. 

Instead, we can consider models with arbitrary slow-roll parameters in the presence of a non-minimal coupling. 
By performing a sort of rotation in the parameter space, we can extend the results continuously to any value of 
effective v and thus to any value of slow-roll parameter e. Recall the scaling (JTTJ) , which relates the that determines 
the background expansion, to the parameter v that governs the behavior of the scalar field mode functions in (|27p . 
Using this, one readily deduces that for a given i/, the slow-roll parameter will have the value e at asymptotically late 
times if the coupling is set to 

? 24(2 -e) ' ( ' 

This means also that we have an asymptotic de Sitter expansion for any v when £ = (9 — 4// 2 )/48. For v = | this 
already requires a rather large coupling £ = — 1/3. It seems thus most natural to consider quasi-de Sitter spaces 
expanding about v = 3/2, since then we need only an infinitesimal coupling £ to have e deviating infinitesimally from 

zero:£=ie-iEr= 2 (|) n . 



This seems to suggest that the result ([45]) is robust to small deviations from exactly exponential expansion. A 
nonzero e, achieved by just adjusting £ e/4, nothing but slightly changes the coefficients of the decaying terms in 
the propagator, whereas the leading logarithm behavior is unaffected. If e is constant, the expansion is described by 

- (1 - e)U = - . (48) 
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FIG. 4. Evolution of the effective equation of state Wbr of the backreaction energy density as a function of conformal time in 
units of | -y | . The backgrounds of these models were displayed in figure[T] Now \~ffc\ = \~f 2 \ = fi. Left: The case of real 7, where 
K — 0. Right: The case of imaginary 7, corresponding to negative curvature K — 2j 2 . In both cases, the backreaction density 
scales as radiation in the early time. The effective equation of state typically diverges, as the backreaction energy density 
changes its sign from positive to negative. Then the scaling approaches a constant again, which is Wb r — > —1 when v — 3/2 
and Wbr —j- —1/3 for the other examples. 



Plugging this into (|37p . we obtain 
//-^A(x; x) 



7*(l-<0 
Ana 



2 r 



D-A 



2 log 



2„2 



4tt/z a 



| 7 2 -(i-k)^l 



(49) 



where £ is given by (|47l) . If e is positive, the prefactor decays, and in super-inflating case, e < 0, it grows. 



IV. BACKREACTION 



In this section we first compute the general expression for the one- loop contribution to the trace in IIV Al This is 
used to analyze the backreaction in inflationary cosmologies. In HVBI we compute the asymptotic contribution to the 
effective energy density. In IIV CI we consider models with contracting and bouncing geometries. 



A. One-loop contribution to the trace 



To compute the expectation value of the trace ©, we have to act with the box on the regular part of (13"T|) and add 
the contribution from the regularization (|43|) . Finally, we must collect the remaining terms from (|34[) . The result is, 
setting D = 4, 



(n\T\n) = - 



(1 - 4v 2 ) (3*1 + 2(1 + h)Hr]) + (2H 2 - (1 + ki)M) + 4 (2H 2 - (1 + k 2 )H') t|V 



327r 2 aV 

[(1 - 4u 2 ) (3 + 2Hr/ - H'rj 2 ) - AH' iW] log a 
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87r 2 a 4 ^ (n- l)r(i/+ \ -n)( 7K ??) 2 " v V 2 + V 



f„r(i/+i + n) /2n+l m W 



—7- \2H 4 + 3H 2 H' + 2{U') 2 - K(4H 2 + m' - AK)\ + —3- M" - 6H 2 W - 2(H') 2 + 2H'K) . (50) 
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We have used shorthand notation for the dimensionless constants k± and k%. In flat and open universes they are 



= -1+7E + log (J*?) -2(1-5C) 
In closed universes the appropriate constant are 



when K < . (51) 



fci = -2 + 7 E + log(^) + S-i#x 2 F 2 (|,2;3,3;-i 
fc 2 = -| + 7£ + log (i^J - 2(3 - 5£)% + SI x 3^2 (f , 1, 1; 2, 3; 



IK 



when If > . (52) 



We also introduced the coefficients f n , denned as 



f 1 when K < , 

r " = ^ ^^x^fn-Litnini-S) when if > . (53) 



IK 



Finally, in the last line of (|50[) the contribution from counterterms is given as 

a mf = 12(D - 4)a 2 , (54) 
a fm = 24(3a{ m - a{ in ) , (55) 

in terms of the coefficients specified in the lagrangian (14TJ1) . As expected, the result ((511)) diverges when 7^ — > 0. The 
divergence is logarithmic for v = 3/2, and power-law ^7^- " for v > |. 



We can obtain the backreaction energy density by integrating (|50|) . and then immediately deduce also the back- 
reaction pressure. This is because we know that the effective quantum fluid must obey a continuity equation similar 
to (|A9[) . and furthermore the trace is T — —pbr + 3pb r , the lower index br referring to the backreaction. Hence we 
have that 

(a 4 Pbr )' = -a 4 HT. (56) 

The integration constant corresponds to a freedom to add a radiation-like component. In the following we show 
numerical results for the integration. In all the numerical examples presented here, we use the minimal subtraction 
scheme and set = 0. Unless otherwise specified, £ = 0, I7I = |7k| = [i, v = 3/2 and K = 0. In realistic models, 
there is a large hierarchy between K and the scale of inflation, Ho in Eq. (|4"4"]) . In units of the Planck scale, a 
GUT-scale inflation for example corresponds to \/8ttGnHo ~ 10~ 5 , where Gm is the Newton's constant, whereas the 
present observational bounds on the curvature radius of the universe, quoted in the first paragraph of this paper, 
translate into SnG n\K \ < 10 -43 . It is reasonable to assume that the parameter j 2 and thus also 7^- are roughly 
on the order of K . The scale factor can be arbitrarily rescaled, but we once fkM = Hq at a = 1, this corresponds to 
the choice 

lC2l = W){^o—) > (57) 
in equation (I2TJ1) or (|2"3"|) in section Hi Bl In the example plots here the numerical ratio is \j\/Ho = 10 -35 . 

We describe the relative significance and the behavior of the backreaction energy density by defining the dimen- 
sionless quantities 

8TTGNa 2 pbr Pbr /ccA 

"6r = Tprpi ' w br = • 58) 

in Z Pbr 

The results in terms of flbr are plotted for an asymptotically de Sitter model in Figure [5] The effective equation of 
state is plotted in the figure |H The background expansions in this figure are precisely the ones depicted in figure Q] 
We see that for both real and imaginary 7, the evolutions have some qualitative similarities. Namely, initially the 
density scales as radiation (wb r w | ), and in the future will scale tend to a negative constant, which in the case v = | 
is Wbr = — 1- In the other cases, the backreaction scales asymptotically like curvature, uib r — —1/3. Note that the 
formally defined effective equation of state can diverge when the backreaction energy density changes its sign. The 
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FIG. 5. Evolution of the fractional contribution to the energy density, \£lbr\ as defined in (|58[) . from the backreaction as a 
function of the e-folding time in flat models with v = 3/2. Left: In the case of nonminimal coupling, the scaling of the quantum 
energy density is modified. When the coupling is negative, the backreaction grows. Right: Examples with different choices of 
H when K — 2-y 2 . As the renormalization scale appears only in the subleading terms, it does not affect the asymptotic results. 



asymptotic scaling we observe in the plots is not obviously seen from (|50[) even in the minimally coupled case, as 
the would-be-leading order terms systematically cancel at late times. The asymptotically leading corrections will be 
discussed in detail below in IIVBI Nevertheless, it holds qualitatively that the backreaction energy density decays if 
< e < 1 , goes like a logarithm in de Sitter case and grows if the universe is undergoing a phantom expansion with 
e<0. 

The effects of varying the coupling £ or the renormalization scale /x in the case v = | in flat and open models 
are shown in figure [5] With a negative coupling £, the relative backreaction contribution grows as a power-law. We 
will see below that this is a generic feature, in agreement with Ref. [2l|. With a positive coupling, the quantum 
contribution decays. As already seen in figure [SJ in the nonminimally coupled case, flb r tends to a constant. This 
constant is determined below in (|60l) . The renormalization scale affects only the transient dynamics. Typically the 
contribution from £lbr is positive in the early times, but becomes negative at the inflationary epoch. In particular the 
time(s) when the sign changes depends upon the renormalization scale. 



B. Asymptotic behavior at constant deceleration 



Assuming that the background has constant deceleration, it is straightforward to integrate (|50J) exactly to obtain 
the energy density. This assumption is always valid in the far future in the above described modelqj. We use equation 
(|48[) and call Hq = W(_a = 1). To take properly into account a possible nonminimal coupling, we should replace 
v — > u% given by (fi~7)) for the background equations, as discussed in IIII C 2l We obtain then 

Pbr = ^4 (r=%) 4 " J (4 " 2 " 1} [23 + 72kl ~ 8(7 + 6klH + 20l/|] 



8(2^ 5 - 1) [3(4^ 2 - 1)(2^ - 3) - 8(2^ e - l)( 7 rf] log (t3^) ~ 16 ( 2i/ ? ~ l W + 4fc 2 - ^)(lKV) 2 
7 |(2^-1) ( 2H oV Y U( - 2I ^ t (l + 2n)(l-2^ + 2n)I>+| + n) 



64ttV \l-2v i ) ^ ™(n + l)n(n-l)r(j/+i-n)(7K?7) 2 ™ 



5 If taken to be valid at all times, this corresponds to setting 7 = in Eq. dl6l . As noted above in section [III Al then one needs nonzero 
K to regularize the infrared. This would correspond to a fine-tuned situation where the energy density of an effective background fluid 
is 8ttGnPb = K/a? + HH^/a 2 ' 1 , the first term canceling exactly the effect of nonzero curvature to the expansion, and the second term 
driving the constant-deceleration expansion. In anisotropic cosmologies such cancelling results in isotropy of the cosmic microwave 
background radiation |42H and might be explained by a presence of a two-form field |43|1 . 
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(2^-l) 2 / 2H 0V V U( - 2 



r fin 



32tt V V 1 - 2u 6 
a inf ( 1 - i/ £ - 4i/| + 4(1 + 4t^) Jfr/ 2 



11-4^(3^ -2) + 8if?7 2 ) 



^ T (^ 2 ) 2 log^ 

2^ - 1 T)i 

The definitions of the constants ki, fe, r n , a™-' 7 and <5-' m are given in formulas (I5HI55|) . 



(59) 



v, the sum in the third line is truncated already at n = max(2, v — |) 
i term. This can be understood, since this term originates from the 
2^ ^s^v.™-. /i/i «n -r^^^f ^ ^^,-.0^^^+ Thus it 



Note that in the minimally coupled case 
due to the coefficient which kills the n = v 

part of the propagator (|3"4"]l that scales as a~ 2 vi 1 ~ 2u . When Eq. (|4"8"]l holds, this part is in fact a constant 
doesn't contribute to the expectation value (|3"9")l that is given by the derivatives of the propagator. One may expect 
that this would change in the case of a massive field, in which case then Qbr would approach a constant. However, 
here we are restricted to the massless field. In the following we extract the leading order term in ([59]) at small 77 in a 
few special cases. 



Let us first set £ = 0. Then, in the case v = 3/2 the background expansion is asymptotically described by u>b 
The fractional backreaction energy density then becomes 



n 



br 



48tt 



± [-5 + 16( 7 ^) 2 logo 2 + 0(( lK v) 2 )] 



-1. 



(60) 



Thus the quantum effect on the expansion tends to a constant which is completely independent of the model pa- 
rameters and the history leading to the de Sitter expansion. The result that effect is negative and magnitude of the 
constant is proportional to the Hubble rate squared in units of the Planck mass was to be expected. 

Let us then look at the quasi-de Sitter spaces and non-minimally coupled field. Then the above result general- 
izes to 



br 



G N [-(l-e)ff ?/P 



86e(2-e) log a 2 



192tt(1 - e)V 

+ (5 - e) (4 - 9e + 12(1 - e)ek 1 ) + 128 (7 - e(7 - 2e)) a mf - 704(1 - e) 3 a fln + 0{{ lK r)f) 



(61) 



Again the leading-order term in (|61[) is parameter-independent. We now infer that the relative backreaction density 
\flbr\ grows, i.e. Wb r < wb, if £ < 0. In this case the background is super-inflating, e < 0, but the backreaction energy 
density still grows faster than ps- This suggests that the quantum backreaction may hinder the Big Rip that is implied 
for the super-accelerating background, in line with e.g. Ref. [3]. For a positive coupling to curvature £ > 0, we have 
the opposite behavior: e > 0, meaning that the background density decays, but pb r will decay faster. Then \Qb r \ will 
decrease proportionally to a negative power of the scale factor. This is due to the positive non-minimal coupling, 
without which the growth is only log-enhanced. The result (|60"|) confirms that the well-known leading-logarithm 
behavior in quasi-de Sitter spacetimes is robust. The density becomes negative, and thus the effect is to slow down 
inflation, in accordance with, e.g. Ref. [45|,|46j]. One expects that the e -suppression of the effect would not appear at 
higher loops, though in our case this remains to be verified. The above deduced behaviors are illustrated in the figure 
El, where we monitor the evolution further inside the inflationary epoch and plot the logarithm of |f2f, r |- Qualitatively 
these results agree with those in Ref. [2lJ, where the infrared was regulated by matching the fluctuation modes in a 
decelerating and accelerating flat FLRW with constant background equations of state. 



2. v = 



In the minimally coupled case, when v — |, % — > — | which can be modelled by e.g. domain walls. In this case 
we find that 



VL br -> [-576 log a 2 + 9 - 54fci - 4096a m/ + 7045 /m + 0(( 7K ??) 2 )] 



(62) 
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It is obvious that the backreaction energy density dilutes more rapidly than the classical background energy density. 
The term scaling as the pg drops out from (I59[) as explained above, and the dominant term will be the following one, 
enhanced by the logarithmic scaling. We have that wi, r = — | + (3a log a) -1 , according with the figure|H where Wb r 
is seen to approach the scaling of curvature from above. 

In the case of nonminimal coupling, the n — 2 term in the sum appearing in (|59[) becomes nonzero and thus 
dominating at late times. The leading contribution is then 
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N 



br 



■(l-e)H rj\ ^10(l-2e) 



[l + C((7 /f ?y) 2 loga)] whence <0, 



9. 



br 



G N [-(1 - e)H ri}— r? 2 10e(2 - e) 



when K > . 



(63) 



We now deduce as above that if £ < (£ > 0), the quantum backreaction grows (decays) like a power- law. Now 
the magnitude of the backreaction depends on the regularization parameter "fx- However, one deduces that it is still 
always negative for arbitrary curvature. There is no e-suppressed log-enhancement. 



In the case v = | , the background fluid is described asymptotically by wb = — § when £ = 0. Now the slowest-dying 
contribution to (|59p comes from the term n = 2 in the sum. The coefficients have a different form if there is a cut-off, 
according to (|53l) . however similarly as above one sees that now uv — > — 1/3 since 



25# 1 V 
708588tt 2 7^ 
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1 + 6{K) 1-16 4 
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o(v 8 ) ■ 



(64) 



It is clear that the magnitude of this term depends on the infrared regulating term 7^ even when the effect of "/k 
upon the background has been completely washed away. In a closed universe, the effective backreaction density (|64|) 
is always positive, whereas it is proportional to the sign of 7^- in open and flat cases. 



With a nonminimal coupling the leading order term is 
G N [-(l-e) J ff 7?] A 35(2-3e 



27T7^6 



[1 + G(( 7 ^) 2 )] ,when#<0 : 



br 



G N l-(l-e)H 0V }— 210(2 - 3e) 
4 r 
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.when.fi: >0. (65) 



The contribution is analogous to the previous case, and also now turns out to be negative for any curvature K and 
parameter 7. 



4. v = 



When v = |, wb 



Pbr 



Pbr 



■| if £ = 0. Analogously to the above case, we obtain now from the n = 3 term in the sum 
735# 16 77 8 



8589934592tt 2 7^ 
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42949672967r 2 7i 



K < 0; 
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K 



0{r) 10 ), K>0. 



(66) 



Now instead the open and flat models have always positive energy density, whereas the sign depends on the ratio 
'y 2 / K in closed models. The next cases follow the analogous pattern: with increasing v the magnitude of the effect 
gets smaller and decays more rapidly with the expansion. Again, if £ 7^ also the n = 4 contributes and will then 
dominate at late times. Given £ < 0, |f2b r | will increase with the scale factor, but its asymptotic sign is always 
negative. 
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C. Backreaction across bounces 



As the final application, let us briefly consider the nonsingular background evolutions described in section III B 31 
Since these evolutions were constructed in such a way that they lead to inflation, the asymptotic effects at 77 — >• 
will be just as detailed above. However, it is of some interest to monitor the evolution of the backreaction also at the 
contracting phase and across the bounce. 

During contraction W 2 ~ K is a constant since a = cie' 7 ' 7 '. If one winds back towards the past r\ — > —00, 
where the size of the universe approaches infinity, it is easy to show by direct integration of the dominating terms 
in (I5U1) that there the backreaction density scales as p& r ~ a~ 4 (with a magnitude depending on a m f , ci, 7 and K). 
Thus in the asymptotic past the quantum effects are negligible. This also means that, as we approach the bounce, the 
backreaction density will grow. The turning point occurs at some finite 77, and the relative magnitude of the various 
terms in (|50[) depends then on the parameters. In particular, ptr can be both positive or negative. At the turnover 
flbr of course diverges. As ordinary quantum fields can acquire an effectively negative energy densities that inevitably 
become dynamically significant as the universe contracts, no new physics may be needed to obtain a bounce even in 
a flat universe. 

In figure [6] we plot the backreaction energy density in the models that were shown in figure [2] The fluctua- 
tions do not feature any pathologies, as H — > 0. In the left panel of the figure [5] we have negative spatial curvature, 
and in the right panel positive. The evolution is very similar in both cases and confirms the behaviors deduced above. 
In both cases, the backreaction density grows as radiation when the universe is contracting and symmetrically decays 
in the expanding phase. The future asymptotic behavior is as discussed in the previous subsection. It remains to be 
studied in detail whether the violation of the null energy condition at the bounce could be naturally caused by the 
backreaction term that increases during the contracting phase, as the results here seem to suggest. 



V. CONCLUSIONS 



Using the operator formalism, we derived the leading quantum corrections from a non-minimally coupled massless 
scalar field in D-dimensional FLRW background. This well understood framework [3] was for the first time applied 
to realistic cosmological models where the spatial sections need not be exactly flat, and furthermore the universe can 
evolve from one phase to another. The propagator was constructed by integrating over the Fourier mode sum (or 
generalized Fourier modes in curved models). The result (j34|) . which in most cases reduces to the simpler form ([35]) . 
was applied to derive the expectation value of the trace of the renormalized stress energy tensor for the quantum 
fluctuations (f5D|) . Also, the asymptotic energy density associated with these fluctuations (JSHJ) was computed to 
analyze their relevance in particular to inflationary cosmologies. 
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A special care was taken to obtain infrared finite results by physical means instead of automatic subtraction at 
play in dimensional regularization (which we employed to get ultraviolet finite answers). We found that a simple 
way to obtain meaningful results is to allow some spatial curvature to exist. A special example of our results is 
that the scalar field in a FLRW universe even with constant deceleration becomes regular in the presence of any 
nonzero curvature (this may also occur for large enough nonminimal coupling, see Figure 1. of Ref. |21|). In general, 
both positive and negative spatial curvature can regulate the infrared. A negative K modifies the evolution of the 
fluctuation modes, in a similar way as a fluid responsible for an early decelerating expansion history, and can thus 
render these modes less singular at very large scales. On the other hand, with a positive K the very large scale 
modes are absent, because the eigenmodes of fluctuations are discretized in a closed universe and the perturbation 
wavelengths cannot be arbitrarily large. 

A parameterization was introduced which allows to study analytically a wide variety models with contracting, 
expanding or bouncing scale factors evolutions. In particular, the parameterization describes realistic inflationary 
expansion histories where inflation was reached from a contracting, decelerating or curvature-dominated era. It was 
found also that by taking into account such a preceding era eliminates the inflationary infrared divergences. Then 
one does not have to set the modes into a Bunch-Davies vacuum in an accelerating stage in the asymptotic past, from 
where any physical trajectory could not have reached the observer, due to the fact that past-eternally inflating space- 
times are not geodesically complete [3l|. Concerning specifically bouncing scenarios, our analysis shows explicitly 
that by assuming a bouncing mechanism, one may not avoid only the initial singularity and geodesic incompleteness 
associated with the background but also the infinities appearing in the fluctuation spectra (we also found that the 
regular backreaction density can grow and have a negative sign during the contracting phase, thus perhaps providing 
the mechanism for a bounce). 

In this set-up then the results have straightforward physical interpretation and furthermore are not affected by 
any ad-hoc regularization that has to be invoked in too simplified cosmological models. We computed the leading 
order backreaction focusing on inflationary cosmologies. In the minimally coupled models, the backreaction tends 
to scale like radiation in an early expanding phase, and then scale like curvature during the inflationary phase. Ex- 
ceptions are the asymptotic de Sitter case, where the backreaction is enhanced logarithmically with the scale factor, 
and super-inflation, where the backreaction can be power-law enhanced. The details of rcnormalization affect only 
the transient dynamics. There sign of the induced backreaction density depends also on the background parameters, 
in particular the spatial curvature. However, in the cases we found the increasing backreaction during inflation, its 
effect was always to slow down the expansion. In the non-minimally coupled case, a generic result is that a negative 
coupling parameter sets the backreaction density growing (relative to the classical background density). The main 
effect of non-minimal coupling is that the field " sees" the background expansion differently, through ([T7]) , otherwise 
the coupling modifies only some finite constants in the results. 
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Appendix A: Field equations for D-dimensional FLRW metric 



The line element on a curved D — 1-dimensional sphere can be written as 

(KId-x = hijd^d^ , (Al) 
where ipi € [0, n] when i — 1, . . . ,D — 2, and i^d-i € [0, 2ir), and the metric is given by 

fkj = -^diag (l,sin 2 (-0i),sin 2 (-0i)sin 2 (V> 2 ), . . . , sin 2 (-0i) sin 2 (V> 2 ) . . . sin 2 (ipn -2)) ■ ( A2 ) 

The role of curvature becomes more obvious when we transform the first angular coordinate into the radial coordinate 
r. These are related by sin(^i) = \f~Kr. Using this relation, we note that the metric in the coordinates (JXJ) is 
generalized to D dimensions straightforwardly as 

g^dx^dx" = a 2 (77) (-dif + 1 + r 2 dn D ^ , (A3) 
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with fi_D_2 given by (|A1|) . In order to determine the Friedmann equations in the curved FLRW geometry, as usually 
we first calculate the nonvanishing Levi-Civita connection coefficients, 



r°oa=H, r;, U5), r°- = -%sin 2 (^)...sin 2 (^-i) 



H 

T 

r i =r * _ Jcot(V>j) if J < « — 1 r j _ f -cot(V>i)sin 2 (^)...sin 2 (V>i-i) ifj<i-l , 

L ij- L ji-\ o if j > i ' » \ if j > i ' { 1 

No sum is indicated over repeated indices above. These coefficients allow us to determine the elements of the Ricci 
tensor and the Ricci scalar 

a 2 R m = {D — i)H'gm ; = , a 2 R l3 = [(£> - 2) ("H 2 + K) + W] gij . (A5) 

a 2 i? = (£> - 1) [(£> - 2) (H 2 +K)+ 2H'] . (A6) 
The elements of the Einstein tensor follow immediately, 

a 2 Goo = ~\{D - 1)(D - 2) {% 2 + X) 5o0 , G« = , 

a 2 G y - = -i(D - 2) [(£> - 3) (H 2 + Jf) + W] gij • (A7) 

We shall assume that the background stress energy tensor is an ideal fluid and thus can be written as T^ v — u^u v (pb + 
Pb) + g^PB, which in the fluid rest frame, specified by u M = (—a, 0, 0, . . . , 0) reads just 

t oq = -PBgoo , Tg = , Ty = -p B gij ■ (A8) 

The divergence of the stress energy tensor vanishes, implying the continuity equation 

p' B + (D-l){ PB + PB )H = 0. (A9) 

It is now straightforward to write the background field equations. The (00) and (ij) components imply the Friedmann 
equations 

where Gn is the Newton's constant. The condition for acceleration is that H' is positive, where we can easily solve 
from the above pair of equations that 

_ 8nG N a 2 [(D - 3)p B + ( D - l)p B ] 

(£>-l)(D-2) ' 1 ' 

From this we see that, just as in four spacetime dimensions, the Universe's acceleration is independent of the curvature. 
Furthermore, to obtain acceleration we require a sufficiently negative equation of state for the background fluid 
wb = Pb/ Pb < — (D — 2>)/{D — 1), which violates the strong energy condition that stipulates w > — 1/(1? — 1). Next 
we use Eqs. (|A10p and (| Al 1|) to rewrite Eq. (|A6|) in terms of the background fluid as 

R=^f(l-(D-l)w B ) PB . (A12) 

Finally we note that though we assumed a positive spatial curvature in (|Alj) . the formulae (|A5|) - (|A12|) are valid for 
arbitrary K. 



Appendix B: Eigenfunctions of the Laplacian on a (D — l)-sphere 

The normal modes of fluctuations in FLRW spaces were discussed in detail in e.g. [13]. Here we sketch the derivation 
in higher dimensions in order to make our conventions explicit and to clarify the interpretation of the wavevector k 
through which we sum/integrate to obtain our real-space propagators. Anticipating possible future extensions of the 
computations at hand to higher loops, we also derive explicitly the radial parts of the eigenfunctions. These will be 
needed in evaluating the propagators at off-coincidence. The angular parts are independent of K (and their square 
in the mode sum may be integrated over using 8.411.7 in [38j). 
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1. Closed universe 



We consider the D — 1-dimensional sphere (IA1|) with the metric given by (|A2|) . We have first assumed that the 
spatial curvature is positive. The Laplacian operator derived from this metric is 



D-l 



s D 



1 = K L „ ; „J-l-n 



sm 



(Bl) 



It is known that on a flat sphere the eigenfunctions of the Laplacian are (higher dimensional generalizations of) 
spherical harmonics, 

Vi^yf^i, . . . , Vr>-i) = + D - 2)^(^1, ■ ■ •, Vd-i) ■ (B2) 
Here M is a collective index for the multipole modes corresponding to D — 3 polar and one azimuthal angle. 

We want to find the eigenfunctions Q(k, x) of the Laplacian (|B1|) . i.e. the solutions to the Helmholz equation 

(V| D _i + fc 2 )Q(k,x) = 0. (B3) 

Since merely the radial dependence changes with nonzero curvature in the metric (|A3|) . the general form of the 
solutions will be 



Q(k,x) = /|(r)^ M (V2,...,V'D-i), 
where the label k will depend on the curvature. The convenient choice will turn out to be 

2 



k 2 =k 2 



D-2 



K. 



The Laplacian decomposes into 



V 2 

v S D 



(B4) 



(B5) 



(B6) 



We are then ready to write down the equation for the radial part of the eigenfunction. In terms of the ip\ , it becomes 



1 n ■ D-2, , x k 2 £(£ + D-3) 

■ d-2,, : d ^ sm m)9^ + -^ + — . 2 / ; \ 

sm (^l) K sm (Vi) 



/^i) = 0. 



The solutions to this equation can be written in terms of the associated Legendre functions, 

/ fc £ W>i) = sin-^(^) [ C1 P£ (costy-i)) + c 2 Q» (costy-i))] , 

where 

1 



A = 



D- 2 



A: 



K 



1 



H = *+-(£> -3) 



(B7) 

(B8) 

(B9) 
(BIO) 



At least at D = 4, it occurs that is a half integer, and we can use, instead of P£ and the pair P x M as independent 
solutions, see formula see 8.737.1 in [38]. By requiring regularity at the origin, using formula 8.756.1 of (38|, we are 



left with only P~^ . Thus 



fM = Nlsm- — (i> 1 )P^ (cos(V^i)) . 



where N[ is a normalization constant. The normalization may be fixed as 



N . 



(2fc + i)L(fc/yx -/i + ±) 

\| 2T(k/VK + [i+ i) 



(Bll) 



(B12) 
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corresponding to the measure 



-l 



dxsm D - 3 (x)f k {x)f k ,(x) = 6 k , k ,5w . (B13) 



Periodicity requires that k/yK is an integer in D = 4, and its lowest eigenmode is k = Dy/K/2. This can be seen as 
follows. The eigenfunction /|(^>i) is single-valued when (see 8.737.2 in |38j). 

f k = cos [(A - /zH/fcW'i) • (B14) 

It follows that k 2 1 K = n(n + 1), where n is a positive integer. The fundamental mode is n = 1, and the values of k 
and £ are then discretized as 

-^L =n+ -(D-2), -n<£<n, n = l,2,.... (B15) 
K 2 



A function may then be expressed as 

oo ^ 

F(rj, . . . , Vd-i) = E E F fo ( n + o( D - 2 ))^ )/L«=a )v ^(^i)«M(^)yf(^2, • ■ • , V>d-i) ■ (B16) 



n=l LM 



Replacing the sum with by integral is an approximation that becomes increasingly good in the UV. We may then 
write the Fourier transformation (a bit schematically) as 



1 f°° 

F(rj, x) = / d D - l kF(rj, k)Q(k, x) 

(2tt)— J§ Vk 



(B17) 



In terms of the wavemode k, shifted with respect to the eigenmode appearing in (|B2|) . the full d'Alembertian acting 
on a scalar function F gives yields 

OF = ( F" + (D — 2)UF' + k 2 F - -(D - 2) 2 KF ) , (B18) 
a z \ 4 / 

leading to our result ([7]) in the case of non-minimally coupled scalar field. 

2. Open universe, K < 

With negative curvature, the analysis is quite analogous. The results can be read from above by analytic continu- 
ation. In particular, the mode functions become then 



and thus 



ftm = J (2fc + * D -2)-i) sinh _^ s)-. 

y dx smh D - 3 (x)f k (x)f k , (x) = S(k - k')St,i, . (B20) 

Now the spectrum is continuous as there are no periodicity conditions. Now k > (D — 2)y/—K/2, thus k > C0. 
The upshot is that by using the index k, we can consider all positive modes as relevant and don't need to introduce 
a nonzero IR cut-off. When matching a particular wave-mode to a physical observable, the subtle question would 
arise, which is the effective physical wavelength of a perturbation, 27r/fc or 2tt jkl However, in the present study this 
ambiguity does not arise as we are considering real-space backreaction quantities, integrated over all Fourier modes. 



Though the basis is complete for these values, and larger modes are not generated by a scalar field, it is interesting to note that the 
most general homogeneous Gaussian random field involves the modes up from k > 0. This has been clarified in [48| . 
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Appendix C: The propagator when K < or j' 2 > 

Here we consider the propagator, barring the case that K > > 7^-. Equation (|35|) can then be written as 

zA(z;a;) = A u -6{K)i8A v , (CI) 

where 

• a 7g' 2 ^ r(l-f + n)I>+l+n) 1 

" a^'V? ^ n!I>+±-n) ( 7 ^) 2 " ' 1 j 

and 5A W appears in closed universes K > to subtract the supercurvature modes. The general expression for it can 
be written in terms of the incomplete beta functions 

•AA (-l)^ 7 g- 2 ^ F(l + n)I> +_§+n) /l 1 \ 1 

In the following we evaluate exactly the propagator at a few special values of v. 

= , ,i ( !_~nt-n, (7l^ 2 -D + 2), (C4) 



iA l = ,0 J ( 2 d'tL 4 W - 3(D - 2) 7 |r7 2 + 3(D - 4)(D - 2)) , (C5) 
2 (2a) u 7r 7k- »y 



zA i = ,0^-2 5 V-z? 6 W - 6 ( D - + 15(D - i)(D - 2hW - 15(D - 6)(D - 4)(D - 2)) . (C6) 



r(i-f) 

(2a) B - 2 7r B 7^T u ?/ 6 

Then consider the terms which take into account the cut-off. In D = 4, Eq. (|C3[) becomes as follows 
•r a Ki ^j r(l + n)r(,+ l + n) / 3 1 5. AK\ 1 

^ " 64^ 7iffl 4 nlr („+!_„) 2 ^ \2' 2 + n ' 2' j*. J ( 7 i^) 2 " ' 1 J 

In the case 1/ = | this reads 



(( 7 | + 4K) r? 2 - 2) - 3 ( 7 | + 4Jf) ((jkv) 2 ~ 2) sinh" 1 (2^ 

iSAs = O* (C8) 

1024 ( 7 | + 4iT)7riaV 

The next couple of cases are 

6^ ((<& + 4^ 7 k) 2 V 4 - 6 7 fc ( 7 | + 4tf) r? 2 + 24X) - 3 7 | ( 7 | + 4if) ^£ + V ( 7 | ^ - 6) sinh" 1 (^) 

10247^(7^+4^)^ + 1^(^)4 

(C9) 



z<5A7 



1024 ( 7 | + 4Kj K ) TT^a^T] 6 



6VK^ + l(rt ( 7 | + 4if) V - 12 7 ^ ( 7 | + 4if) V 



120 7 |^ (t£ + 4RT) r? 2 + 120X (5 7 | + 8/-T) J - 3 (7$. + 4iT 7jf ) V {lW ~ 12) sinh^ 1 ( j j (.C10) 



These special cases are analyzed and evaluated numerically in section [TV 



[1] E. Komatsu et al. (WMAP), [Astrophys. J. Suppl. jl80, 330 (2009), |arXiv:0803.0547 [astro-ph]| 



23 



A. H. Guth, |Phys~ Rev. D23, 347 (1981). 

A. A. Starobinsky, |Phys. Letf]B117, 17 5 (1982). 

A. Vilenkin and L. H. Ford, |Phys. Rev.| D26, 1231 (1982). 

L. P. Grishchuk, Sov. Phys. JET P 40, 409 (1975). 

T. Janssen and T. Prokopec, jAnnals Phys] 325, 9 48 (2010), |arX"iv:0807.0447 [gr-qc]| 

N. C. Tsamis and R. P. Woodard, J. Math. Phys. 48, 052306 (2007), arXiv:gr-qc/0608069 



T. Janssen and T. Prokopec, JCAP 0705, 010 (2007), arXiv:gr-qc/0703050 
C. Germani and A. Kehagias, JCAP 0911, 005 (2009), arXiv:0908.0001 [astro-ph.CO] 
T. S. Koivisto and N. J. Nunes, |Phys. Rev.|D80 , 103509 (2009), arXiv:0908.0920 [astr o-ph.CO]] 
N. C. Tsamis and R. P. Woodard 7|Annals PhysTl 253, 1 (1997), [arXiv :hep-ph/9602316 

V. F. Mukhanov, L. R. W. Abramo, and R. H. Brandenberger, |Phys. Rev. Lett. 1 7871624 (1997), |arXiv:gr-qc/9609026 
T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. Lond. A357, 381 (1977). 
N. Birrel and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge, UK: University Press, 1982). 

A. Vilenkin, |Nucl. Phys] B226, 5 27 (1983). 

N. C. Tsami s and R. P. W oodard, |Class. Quant. GravTI ll, 2969 (1994). 

A. D. Linde, |Phys. Lett] B 114, 431 (19 82). 

B. Losic and W. G. Unruh, [Phys. Rev.|D74 023511 (2006), [axXiv:gr-qc/0604T22| 

B. Losic and W. G. Unruh, [Phys. Rev. Lett] 101, 111101 (2008), |arXiv:0804.4296 [gr-qc]| 

T. Janssen and T. Prokopec, Class. Q uant. Gray. 25, 055007 (2008)7 |arXivT07 07.3919 [gr-qc] 

T. M. Janssen and T. Prokopec(2009), |arXiv:0906.0666 [gr-qc] | ~ 

T. M. Janssen, S. P. Miao, T. Prokopec, and R. P. Woodard, |Class. Quant. Grav.| 25, 245013 (2008), 
|arXiv:0808 .2449 [gr-qc] | 
A. D. Linde,|Phys. Lett]B116, 335 (1982). 



S. B. Giddings and M. S. Sloth(2010), arXiv: 1005. 1056 [hep-th], 

C. T. Byrnes, M. Gerstenlauer, A. He becker, S. Nurmi, and G. T asinato(2010), |arXiv:1005.3307 [hep-th]| 
L. Senatore and M. Zaldarriaga(2009), [axXiv:0912.2734 [hep-th]| 

E. O. Ka hya, V. K. Onemli, a nd R. P. Woodard (2010),|arXiv: 1006.3999 Tastro-ph.CO]| 

D. Seery, pass. Quant GravT] 27, 124005 (2010), |arXiv:1005. 1649 [astro-ph-COF " 

N. Arkani-Hamed, S. Dubovsky, A. Nicolii; E. Trincherini, and G. Villadoro, JHEP 05, 055 (2007), 
|arXiv:0704. 1814 [hep-th] 

M. Novello and S. E. P. Bergliaffa, |Phys. Rept]463, 127 (2008), |arXiv:0802.1634 [astrojph] 



A. Borde, A. H. Guth, and A. Vilenkin, |Phys. Rev. Lett. 1 90, 151301 (20 03), |arXiv:gr-qc/0110012 
T. S. Koivisto, [Phys. Rev.| D82, 04 4022 (2010), |arXiv: 1004.4298 [gr-qc]" 
C. Barragan and G. J. Olmo(2010), arXiv: 1005.4136 [gr-qc] 



G. Calcagni and G. Nardelli(2010), arXiv: 1004.5144 [hep-th 
T. Biswas, T. Koivisto, and A. Mazumdar(2010), |arXi v: 1005.0590 [hep-th] 
T. Biswas, A. Mazumdar, and A. Shaneloo(2010), |arXiv:'l 003.3206 [hep -th] 
T. Biswas and A. Mazumdar, |Phys. Rev.] D80, 023519 (2009), ]ajXiv:0901.4930 [hep-th][ 

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products (New York, USA: Academic Press, 1965) 

S. P. Miao, N. C. Tsamis, and R. P Woodard(2010), |arXiv: 1002 4037 [gr-qc]| 

V. K. Onemli and R. P. Woodard, |Class. Quant. Grav.|19, 4 607 (2002)" |arXiv:gr-qc/0204065| 

E. O. Kahya, V. K. Onemli, and R. P. Woodard, |Phys. Rev] D81, 0235 08 (2010), |arXiv:0904.4811 [gr-qc]1 
S. Carneiro and G. A. Mena Marugan, Phys. Rev. D64, 083502 (2001), arXiv:gr-qc/0109039 
T. S. Koivisto, D. F. Mota, M. Quartin, an d T. G. Zlosn ik(2010), |arXivTl0 06.3321 [astro-ph.CO] 
S. Nojiri, S. D. Odintsov, and S. Tsujikawa, |Phys. Rev.| D71, 063 004 (2005), |arXiv:hep-th/0501025| 
N. C. Tsamis and R. P. Woodard, |Nucl Phys.| 'B474, 235 (199 6), |arXiv:hep -ph/9602315 



L. R. W. Abramo, R. H. Brandenberger, and V. F. Mukhanov, Phys. Rev. D56, 3248 (1997), arXiv:gr-qc/9704037 
E. R. Harrison, |Rev. Mod. Phys.|39, 862 (1967). 



D. H. Lyth and A. Woszczyna, Phys. Rev. D52, 3338 (1995), arXiv:astro-ph/9501044 



